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FOREWORD 


The  Energetic  Materials  Division  has  an  ongoing  program 
investigating  double-base  propellant  combustion  mechanisms.  This 
effort  is  supported  by  the  NAVSEA  contract  N000248WR10192 
element  6115 3N  and  by  NSWC  internal  Independent  Research  funding, 
1R01AA230. 

In  support  of  these  investigations  the  problem  arose  of 
characterizing  two  closely  spaced  chromatographic  peaks.  This  paper 
reports  the  use  of  the  Scofield  unfolding  algorithm  as  a  method 
of  mathematical  deconvolution.  Due  to  the  Scofield  algorithms 
property  of  preserving  the  sign  of  intermediate  results  a  superior 
end  result  is  achieved.  The  Scofield  algorithm  was  initially 
developed  at  the  U.S.  Naval  Radiological  Defense  Laboratory  and 
was  further  applied  here  at  the  Naval  Surface  Weapons  Center  to 
various  radiological  programs.  This  report  represents  a  further 
application  of  this  remarkable  algorithm.  No  endorsement  of 
commercial  products  mentioned  in  this  report  is  implied. 

The  authors  express  their  appreciation  to  Cindy  McPherson  for 
assistance  in  preparation  of  the  graphics,  and  to  representatives 
of  Waters  Associate  and  Perkin-Elmer  Corporation  for  helpful 
comments. 
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INTRODUCTION 


Satisfactory  resolution  of  chromatographic  peaks  is  necessary 
for  both  qualitative  and  quantitative  analysis.  Despite  the  wide 
variety  of  possible  combinations  of  detectors,  columns  and  solvents, 
it  is  common  to  encounter  poorly  resolved  or  unresolved  peaks. 

Various  methods  of  mathematical  analysis  of  the  curves  can  then  be 
applied  as  a  last  resort. 

One  such  method  models  the  chromatogram  as  a  summation  of 
gaussian  curves  or  modified  gaussian  curves  by  least  squares 
techniques^ .  This  method  can  be  useful  but  suffers  to  the  extent 
that  it  does  not  handle  non-gaussian  or  tailing  peaks.  Another  broad 
class  of  techniques  is  digital  filtering  in  one  or  more  of  its 
disguises.  For  example,  if  the  shape  of  an  unblended  peak  is  shown 
it  is  assumed  to  represent  the  instrumental  response  function.  An 
"optimal"  filter  can  then  be  accomplished  by  dividing  the  Fourier 
transform  of  the  unresolved  data  by  the  Fourier  Transform  of  the 
response  function2.  The  advent  of  the  Fast  Fourier  Transform  makes 
this  elegant  technique  efficient  and  readily  applied  to  automated 
chromatographs.  However  the  optimality  of  the  Fourier  filtering 
process  depends  upon  the  fundamental  assumption  of  mathematical 
linearity.  Although  the  physical  peaks  do  add  linearly,  it  is  also 
true  that  a  negative  amount  of  matter  is  a  physical  impossibility 
(i.e.  a  non-linear  constraint).  Thus  we  have  the  a  priori  knowledge 
that  the  solution  to  the  resolution  problem  can  only  have  positive 
peaks  when  the  specimen  affects  the  detector  positively  and  only 
negative  peaks  when  the  specimen  affects  the  detector  negatively.  The 
Fourier  Transform  admits  the  possibility  of  solutions  of  the 
opposite  sign,  consequently  optimality  can  not  be  proven. 


^Scott,  Charles  D.,  Chilcote,  Dennis  D.,  Pitt,  W.  Wilson,  "Method 
for  Resolving  and  Measuring  Overlapping  Chromatographic  Peaks  by 
use  of  an  On-Line  Computer  with  Limited  Storage  Capacity",  Clinical 
Chemistry  16  (8)  637-42  (1970). 

^Rabiner,  Lawrence  R.  and  Gold,  Bernard,  Theory  and  Application  of 
Digital  Siqnal  Processing,  Prentics  Hall,  Inc.,  Englewood  Cliffs, 

NJ  (1975)  . 
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Although  no  claims  to  optimally  are  known,  the  iterative 
unfolding  technique  of  N.  E.  Scofield^  admits  only  positive  solutions, 
is  easy  to  implement,  and,  as  will  be  shown  here,  gives  excellent 
results.  The  Scofield  algorithm  was  designed  to  solve  the  more 
complicated  problem  of  two-dimensional  unfolding  while  we  require 
here  one  dimensional  unfolding  otherwise  known  as  deconvolution. 

The  two  dimensional  problem  occurs  when  the  response  function 
itself  it  not  constant.  While  this  is  true  of  chromatograms  in  the 
most  general  case  we  limit  ourselves  here  to  deconvolving  only  a 
portion  of  a  chromatogram. 

This  allows  us  to  use  the  degenerate  one-dimensional  case  of  the 
Scofield  algorithm  which  has  the  additional  advantages  that  it  can 
be  programmed  into  any  microcomputer  which  possesses  the  following 
capabilities:  (1)  ability  to  store  digitized  detector  output  data  as 
a  function  of  time  in  an  array  of  a  length  equal  to  the  number  of 
equidistant  time-slices  taken,  (2)  minimal  BASIC  logic  (IF/THEN, 
FOR/NEXT,  GOTO,  LET,  +,  -,  /,  *)  or  equivalent,  and  (3)  enough  RAM 
space  to  contain  four  detector  output  arrays  and  the  BASIC  program 
itself  (about  50  logical  lines  for  one  no-frills  version  which 
included  data  acquisition) . 

In  our  laboratory  this  algorithm  has  been  programmed  and  run 
on  a  Sigma  115  GC  System  console  (Perkin-Elmer ,  Norwalk,  CT)  and  on 
a  Model  720  HPLC  System  Controller  (Waters  Assoc.,  Milford,  MA) . 
Examples  of  the  BASIC  programs  which  implement  the  algorithm  for 
these  two  instruments  are  included  in  the  Appendix. 


DESCRIPTION  OF  ALGORITHM 


The  following  variables  are  defined  in  the  deconvolution 
algorithm: 

C (x)  =  an  array  which  contains  baseline  corrected  time-slice 
data  for  the  sample  chromatogram  that  is  to  be 
deconvolved. 


R (x)  -  an  array  which  contains  baseline  corrected  time-slice 
data  for  a  reference  chromatogram  of  a  single  pure 
substance,  such  as  one  component  of  the  sample  mixture. 

N  =  the  length  of  C  (x) ;  the  number  of  time-slices  taken  to 
fill  C  (x) . 


N1  =  the  length 
fill  R (x) . 
width  used 


of  R (x) ;  the 
N1  need  not 
for  C (x)  and 


number  of  time 
be  equal  to  N, 
R(x)  should  be 


slices  taken  to 
but  the  slice 
the  same . 


^Scofield,  N.  E. ,  "Iterative  Unfolding",  Applications  of  Computers 
to  Nuclear  and  Radiochemistry,  National  Academy  of  Science  Monograoh 
NAS-NS-3107 ,  G.  D.  O'Kelley,  ed.  (1962). 
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xl  =  an  integer  such  that  1  £  xl  <  Nl,  and  R(xl)  <  R(x)  for  all 
1  <  x  <  Nl  (i.e.,  xl  is  the  number  of  the  time-slice  which 
contains  the  "peak"  of  the  reference  chromatogram) . 

Nl(x)  =  an  array  of  length  N  which  is  filled  by  the  computer  with 
the  time-slice  data  for  the  deconvolved  chromatogram  Nl(x) 
is  initially  set  equal  to  C (x) ,  and  is  updated  with  each 
iteration  of  the  deconvolution  algorithm. 

Cl(x)  =  an  array  of  length  N  which  is  used  by  the  computer  to  store 
intermediate  values  during  each  iteration. 


The  algorithm  defines  the  following  initial  conditions: 

Nl (x)  =  C (x)  for  all  1  <  x  <  N. 

Each  iteration  of  the  deconvolution  algorithm  updates  N(x) 
by  the  following  equations: 

Y2 

Cl(x)  =  Z  R (xl+i) *N (x+i)  for  all  1  <  x  <  N 

(previous  iteration) 


Where:  Y1  =  max  ( (1-x) ,  (1-xl)) 

Y2  =  min  ( (N-x) ,  (Nl-xl) ) 

N(x)  =  (C (x) /Cl (x) )  *  N (x)  for  all  1  <  x  <  N 

(previous  iteration) 


The  vector  N(x)  can  be  renormalized  after  each  iteration  to 
conserve  area  or  peak  height.  The  deconvolved  chromatogram  N(x) 
shows  the  most  marked  improvement  over  C(x)  after  one  or  two  iterations. 
After  this,  further  improvement  occurs  approximately  as  the  logarithm 
of  the  number  of  iterations  performed.  With  each  iteration,  the 
chromatogram  becomes  increasingly  noisy,  so  that  it  is  advantageous 
to  iterate  as  few  times  as  possible  to  obtain  a  sufficiently  sharpened 
chromatogram . 


EXAMPLE 

Figure  1  shows  a  simple  chromatogram  which  was  deconvolved  on 
a  Waters  Model  720  HPLC  System  Controller.  The  sample  was  a  mixture 
of  38%  by  weight  2-nitrodiphenylamine  (2NDPA)  and  62%  diethylphthalate 
(DEP).  The  mixture  was  separated  on  a  Waters  brand  silica  radial 
compression  column.  The  mobile  phase  used  was  3/1,  2,2,4-trimethyl- 
pentane/tetrahydrofuran.  A  254nm  UV  absorbance  detector  was  used. 

The  original  chromatogram  is  shown  in  Figure  1.  The  major  peak 
eluting  at  2.63  minutes  is  2NDPA;  the  shoulder  is  DEP. 
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Figure  2  shows  an  expanded  segment  of  this  same  chromatogram. 

It  is  composed  of  38  data  points  taken  at  0.03  minute  intervals. 

The  deconvolution  process  was  performed  on  these  data.  The  deconvolved 
chromatograms  are  shown  in  Figures  3  and  4  after  undergoing  two  and 
ten  iterations  of  the  algorithm  respectively.  For  convenience,  after 
each  iteration  the  data  were  normalized  so  that  the  maximum  peak  height 
was  exactly  full  scale  on  the  chart.  After  ten  iterations,  the  main 
peaks  are  sharp  enough  that  the  DEP  is  almost  baseline  resolved.  The 
area  of  the  two  peaks  (times  response  factors)  give  an  analysis  for 
the  sample  of  35%  2NDPA,  65%  DEP.  This  is  in  good  agreement  with  the 
sample's  actual  weight  make-up  of  38%  2NDPA,  62%  DEP. 

In  Figure  4  a  shoulder  can  be  seen  on  the  DEP  peak  which  was 
totally  invisible  in  the  non-deconvolved  chromatogram.  With  further 
iterations,  this  shoulder  too  would  become  a  baseline  resolved 
peak.  This  fact  points  out  a  limitation  on  the  power  of  the  iterative 
algorithm:  an  energy  bump  in  the  chromatogram,  whether  signal  peak  or 
just  random  noise,  is  sharpened  with  each  iteration.  Without  further 
prior  knowledge  the  reality  of  such  features  is  an  epistemological 
question  not  addressed  here.  The  main  utility  of  this  deconvolution 
technique  is  in  the  resolution  of  known  peaks.  It  is  at  best  a  very 
risky  tool  for  detecting  previously  unseen  peaks.  One  can  however 
construe  such  low  level  peaks  as  a  hint  directing  future  work. 


10 


NSWC  TP.  81-307 


0  PER -.TOR  :  DORN  CARLSON 

SAMi  i.E  SAME :  SAMPLE  <  2-NDPA  PLUS  DIETHYL  PHTHALA'S''  ACTUAL  L1'  TRACE 

SAMPLE  PC’S  01/  I i  1 J  -UOLUME  0003/  NO  OF  INJ  i  /  RUN  TIME  00=04/  SM 
HISP  CODES  GENERA  TED =  Sj-SIUL 


MAY.  26/  1381  13=31:25  CHART  1.00  CM/MIN  FLCW  _  J 

SUN  #3  _ 

rni  iimn  SOL  LENT  UPR  IJ: 


2.0«  ML/M  IN 


COLUMN  bULVtrn 

EXTERNAL  STANDARD  QUANTITATION 

PPAK#  AMOUNT  RT  EXP  RT 

36604.10000  2.83  36®04*S3  L 

TOTAL  36604.19000 

PEAK#  AMOUNT  RT  ’• 

36604.10000  2.63  100.00 

‘  OTf.  ,■  3660 4.10000  100.00 


RF 

0 . 000000E0 


Figure  1  Liquid  Chromatogram  of  2NDPA/DEP  Mixture 
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Time,  Minutes 


Figure  3  Deconvolved  Chromatogram  (two  iterations) 
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APPENDIX  A  -  "PROGRAM  DOCUMENTATION" 


The  following  BASIC  programs  are  deconvolution  programs 
written  to  operate  on  the  Waters  Associates  programmable  integrator 
the  Model  720/730,  and  on  the  Perkin-Elmer  Sigma  115  instrument. 
Figures  A-l  and  A-2  are  notated  excerpts  of  the  programs  which  are 
specifically  concerned  with  the  deconvolution  iteration.  Figures 
A-3  and  A-4  are  complete  programs  which  include  some  aspects  of 
data  acquisition,  baseline  corrections,  and  replotting. 
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7080 

7100 

7102 

7105 

W10  8 

'*7120 
^7140 
>7160 
57180 
w7l85 
>7190 
w  7260 
o7280 
£7290 
m7300 
w  7320 
7340 
7360 
7380 
7400 
7410 
7420 
7460 
7500 
7560 
7570 
7580 
7600 
7620 
7640 
7660 
7680 
7700 
7720 
7740 
7760 
7780 
7800 
7820 
7840 


CPRINT  “DECGNUOLUT ING. . . .  " 

CPRINT 

FOR  X=1  TO  N  '  Set  up  initial  conditions,  N(x)  *  C(x) 

N<  X  )=C(  X  > 

NEXT  X 
CPRINT 
N0=0 

D-CHR( 13 > 

CPRINT  “NO.  OF  ITERATIONS  =  ",N0  , 

D=CHR< 7  ) 

'  Do  a  total  of  ten  iterations 


'  Deconvolve  each  point  on  curve,  x 
*  from  x  =  1  to  x  =  N 

1  Set  Cl(x),  intermediate  array,  *  to  0 
'  Y1  =  max(  1-x,  1-X1  ) 

Y2=N-X 

IF  <  Nl-Xl )<<  N-X  )  THEN  Y2=N1-X1  •  Y2  =  min(  N-x,  Nl-Xl  ) 

I=Y1-1 

1=1+1  '  Do  summation  from  i = Y1  to  i = Y2 

Cl(  X )=C1<  X  >+R'  XI  + 1 >*N<  X+I  ) 

IF  I< Y2  THEN  7410  '  Next  i 

IF  X<N  THEN  7290  '  Next  point  on  curve  (next  x) 

X=0 
Ml  =0 
X=X+1 

IF  C1<X>=0  THEN  7660  '  Division  by  0  excluded 

N(  X  )=N<  X  )*C<  X  )/Cl<  X  ) 

GOTO  7680 
N<  X )=0 

IF  N< X )>M1  THEN  M1=N<X)  '  Find  Ml,  the  maximum  height  of  curve 

IF  X<N  THEN  7580 

X=0 

X=X+1 

N<X)=INT(N<X)*100000/Ml+.5)/1000  '  Divide  each  point  by  Ml/100 

IF  X<N  THEN  7740  '  (height  normalization) 

N0=N0+1 

IF  N0=2  THEN  7860 
GOTO  7160 


IF  N0S10  THEN  8000 
REM  "OECONOOLUTION" 

X  =  0 
X=X+1 
Cl<  X  )=0 
Yl=l-X 

IF  X1<X  THEN  Y 1 = 1 — X 1 


Figure  A-l  BASIC  Subroutine  Used  to  Deconvolve  Chromatograms  on 
Waters  Model  720 


A- 2 


NSWC  TR  81-307 


3 o o  LET  S8=® 

7-m  i  c-r  m_a 
ji  <i  uC  <  i  i—V 

333  FOR  X=1  TO  N 

398  IF  Jt<N<X>  THEN  LET  !1=N<X> 

438  NEXT  X 

413  FOR  X=1  TO  N 

423  LET  N<X>=N<X>*78^« 

433  NEXT  X 

463  FOR  X=1  TO  INT<N^2> 

473  LET  81  =  1 NT <N<2*X-1 ) ) 

433  LET  R2=INT<N<2*X>> 

IF  81=82  THEN  PRINT  T88<R1>I *i 
IF  81 <82  THEN  PRINT  TR8<81>;“' 
IF  81)82  THEN  PRINT  T88<R2>;**, 


493 
533 
o  1 3 


Find  M,  maximum  height  of  curve 


'  Divide  each  point  on  curve 
'  by  M/70  to  normalize 

1  Lines  460-530  trace  out  curve 


;T8B<8i>; 


c-.rt  w 

itCn  •  A 

'2)INT<N/2) 


Do  ten  iterations 


533  IF  N>'2;INT(N/2)  THEN  PRINT  TRB<INT<N<N> > )  1  “ ' 
543  PRINT  -ITER8TICN  NO.  “1N0 
553  LET  N8=N8+l 
563  IF  N8=l 1  THEN  1830 
573  FOR  X=1  TO  N 
533  LET  D<X>=8 
638  LET  Vl=l-X 

613  IF  X1<X  THEN  LET  Vl=t-Xl 
623  LET  V2=N-X 


633  IF  <Nl-Xl;<(N-X)  THEN  LET  V2=N1 
648  FOR  I=V1  TO  V2 
658  LET  D<X>=D<X>+R<X1+I>*N<X+I> 

663  NEXT  I 
673  NEXT  X 
690  FOR  X=1  TO  N 
738  IF  OQw=0  THEN  738 
713  LET  N<X)=N<X)*C<X)/D<X) 

723  GOTO  748 
733  LET  N(X)=8 
743  NEXT  X 
753  GOTO  373 
1338  REN 


’  D(x)  corresponds  to  Cl(x)  in 
’  the  algorithm 
1  Yl  «  max(  1-x,  1-Xl  ) 

-XI  '  Y2  -  min(  N-x,  Nl-Xl  ) 

'  Do  summation  from  i  -  Yl  to  i  m  Y2 


'  Division  by  0  excluded 


Figure  A-2  BASIC  Subroutine  Used  to  Deconvolve  Chromatograms  on 
Perkin-Elmer  Sigma  115 
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*EM  » DECONU . BAS  DECONUOLUT I ON  5 '6/8  1 " 

:  RE M  “DORN  CARLSON" 

10  DEBUG 
ZO  CLEAR 
38  D=BACKF 
40  D=KILLF 
58  D=PAGEF 
50  D=MODUL£F 
70  D=ERASEF 
4000  DIM  R< 100) 

4020  DIM  0(500) 

4040  DIM  C< 500 ) 

4060  DIM  C 1(500) 

4080  DIM  N< 500  ) 

5000  ERASE 
5010  0=BACKF 
5020  K 1 =0 

5040  CPF  I  NT  "RESPONSE  FUNCTION  (  REFERENCE  PEAK)" 

5060  GOTO  5120 
5080  ERASE 

5100  CPRINT  “RAW  CHROMATOGRAM" 

5120  CPRINT 

5140  CPRINT  "1  =  COLLECT  DATA  AS  SAMPLE  IS  RUN" 

5150  X=0 

5160  CPRINT  "2  =  LOAD  DATA  FROM  TAPE" 

5165  IF  K 1 = 1  THEN  5180 

5170  CPRINT  "3  =  USE  SAME  RESPONSE  FUNCTION" 

5130  C INPUT  ,X 
5200  D=LIN( 7 ) 

5220  IF  X=2  THEN  9000 
5230  IF  X=3  THEN  6100 
'235  CPRINT 

5240  CPRINT  "INITIAL  TIME,  FINAL  TIME,  INTERNAL" 

5230  C INPUT  , TAB( 2 ), R1 , TAB( 16 ), R2, TAB(  28 ), R3 
5230  CPRINT 
5290  CPRINT 

5300  CPRINT  "HIT  (ENTER)  AFTER  RUN  STARTS" 

5320  IF  BACKF=0  THEN  5320 
5340  TO  PPI  ASK  RUN_TIME=CON( Q ) 

5360  Q1=SYSTIM 

5330  Q2=  INK  Ql/1 00  )-40* I NT(  01/1 0000  >+ 1 NT(  <01-1 004 I NT(  0 1  ■■  1 00  )  >/ .  6+ .  5  V 1 0 
0 

5400  R4=02-Q-.O1 

5401  REM  "R4=DIFF.  BTN.  RUN  TIME  &  SYSTIM" 

5410  N=1 

5420  D=CHR( 7 ) 

5430  R=R1 
54‘  i  01 =SYSTIM 

5450  02= I NT( 01/1 00  )-40* I NT( Ql/1 0000  )+ 1 NT( < 01-1 004 1 NK 0 1 / 1 00 ) )/ . 6+ . 5 >/ 


Figure  A-3  Complete  BASIC  Deconvolution  Progrcim  for  Waters  Instrument 
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4 -pi 

5*-  c  Q2\R  THEN  5*40 

550b  3  PPI  ASK  INPUT_LEUEL=CGN<  Q  :• 

5520  .  =y 
5540  .-■=  ■  . 

55 60  P 

553W  *  THEN  56'0 

5600  GOTO  5440 
5620  D=CHR<  7 ) 

5640  OPR I NT  "DATA  COLLECTION  COMPLETE" 

5660  CPRINT  "CORRECTING  FOR  BASELINE  DRIFT 
5630  N=N-1 
5700  A=C( 1 > 

5720  B=C(N> 

5740  X=0 
5760  X=X+1 

5780  C<  X >=C(  X >-< <  X- 1 >*B+< N-X  ) t A >/<  N- 1  > 

5790  IF  C(X><0  THEN  C(X)*0 

5800  IF  X<N  THEN  5760 

5820  IF  Kl  =  l  THEN  6040 

5340  S1=R1 

5360  S2=R2 

5880  S3=R3 

5900  N1=N 

5910  Ml =0 

5920  X — 0 

5940  X=X  +  1 

5960  R<  X >=C<  X > 

5980  IF  RCX)>M1  THEN  X1=X 
6800  IF  R<  X >>M1  THEN  M1=R<X> 

6020  IF  X< N  THEN  5940 
6030  GOTO  6068 
6840  X=0 
6045  Ml  =  l 
6043  X=X+1 

■'050  IF  C<  X  )>M1  THEN  M1=C<X) 

2  IF  X\N  THEN  6048 
^055  X=0 
6056  X*X+1 

->357  CC  X  )=  I NTC  C<  X  >*  1 00000/M  1  + .  5  V  1 000 

.053  IF  X< N  THEN  6056 

6060  GGSUB  6380 

6100  IF  Kl=l  THEN  6160 

- 120  K 1  =  1 

6140  GOTO  5980 


Figure  A-3  (Continued) 
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clod  CPRINT  “RESPONSE  FUNCTION",  TAB<  22  •,  ""RAi'E  QATA" 
6200  CPRINT 


i  •-?  CPRINT  N1 ,  "  POINTS", TAB< 22 ), N, -  POINTS" 

-—40  CPRINT  SI,"  START" , TAB< 22  )■ R1 . "  START" 

-•260  CPRINT  S3,"  INTERNAL" 

'230  CPRINT  XI, "  PEAK  MAX" 

6320  CPRINT 
6340  GOSUB  6330 
6360  GOTO  6440 

6380  CPRINT  "MIT  <ENTER>  TO  CONTINUE" 

6400  IF  BACKF=0  THEN  6400 
6410  ERASE 
6420  RETURN 

6440  CPRINT  “1  =  STORE  RESPONSE  FUNCTION  ON  TAPES  « 
6460  X=0 
6480  C INPUT  X 
6500  CPRINT 


6510  IF  X=0  THEN  7*40 

6520  REM  m**LINES_6"20_TO_7020_NOT  IN 

7040  CPRINT 

7060  GOSUB  6380 

7030  CPRINT  "DEC 


7100  CPRINT 
7182  FOR  X=1  TO 
7105  N<  X  >=C<  X ) 

7108  NEXT  X 
7120  CPRINT 
7140  N0=0 
7160  Q=CHR( 13  ) 

7130  CPRINT  "NO.  OF  ITERATIONS  =  ".NO 
7135  0=CHR<  7  ) 

7190  IF  N0= 10  THEN  8000 
7260  REM  "OECONUOLUTION" 

7280  X=0 
7290  X=X+1 
7300  Cl< X )=8 
7320  Yl=l-X 

7340  IF  X1<X  THEN  Yl-l-Xl 
7360  Y2=N-X 

7380  IF  <  Nl-Xl  )<<  N-X  >  THEN  Y2=N1-X1 
7400  I = Y 1 — 1 
7410  1=1+1 

7420  C 1 <  X  >=C 1 <  X  )+R(  X 1  + 1  ):*N<X+I  > 

7460  IF  I< Y2  THEN  7410 
7500  IF  X< N  THEN  7290 
7560  X=0 
7570  M 1 =0 
7580  X=X+1 

7600  IF  C 1 ( X  )=0  THEN  7660 
7628  N< X )«N< X  >#C< X  )/Cl<  X  > 

Figure  A-3  (Continued) 
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WE Z 


7640  GOTO  7680 
7660  N<X>=0 

7680  IF  N( X  >>M1  THEN  M1*N<X> 

7700  IF  XCN  THEN  7500 
7720  X=0 
7740  X=X+1 

7760  N(X)=INT<N<X)*100000/|11  +  .5)/1000 

7730  IF  XCN  THEN  7740 

7800  N0=N0+1 

7820  IF  N0=2  THEN  7860 

7840  GOTO  7 ‘  '0 

7860  X=0 

7880  X=X+1 

7900  Q< X  )=N<  X  ) 

7920  IF  X<N  THEN  7880 
7940  GOTO  7160 
8000  CPRINT 

8020  CPRINT  "STORING  DATA  IN  GRADIENT  TABLE,  PUMP" 

8040  CPRINT  "SET  13" 

3850  CPRINT 

3060  CPRINT  "RAW  DATA  5£A“ 

3080  CPRINT  "2  ITERATIONS 

3100  CPRINT  "10  ITERATIONS  JjC" 

8120  TO  PUMP  SET  EDIT_SET=CON< 13 > 

3140  TO  PUMP  SET  ERASE_GD_TABLE=1 

3160  TO  PUMP  SET  TOTal_FLOW=CON< 0, C< 1 >,D< 1  ).N< 1  Vi 

3130  R=R1 

8200  N2= 1 

8220  TO  PUMP  SET  GRAOIENT_TABLE=CON<  R  ,  0,  C<  N2  ) .  D<  N2  >  ■  N>:  N2  >  ■  “06" 
8240  R=R+R3 
3268  N2=N2+1 

3280  IF  N2>N  THEN  8320 
3300  GOTO  8220 
3320  D=CHR( 7 ) 

8330  X— 0 

8340  CPRINT  "1  =  STORE  DATA  ON  TAPE?  "  , 

3350  D=8ACKF 

8360  C INPUT  ,X 

8330  IF  X=8  THEN  3580 

3400  CPRINT  "NAME  OF  FILE  <■  , 

3402  D=CHR(  120  )+CHR<  120  >+CHR<  120>+CHR<  120  )+CHR<  120>+CHR<  120  > 
3404  D=CHR(  46  >+CHR<  120)+CHR<  120>+CHR(  120)+CHR<41  >+CHR<32> 

*4?0  C INPUT  -STRING 
40  CPRINT 

50  CPRINT  "< SAME  PUMP  SET  13  ONLY!)" 

30  CPRINT  "HIT  CENTER)  WHEN  TAPE  IS  READY" 

30  IF  BACKF=0  THEN  3500 

Figure  A_3  (Continued) 


A-7 


NSWC  TR  81-307 


.020  TSAUE  PARAMETER,  STRING 
3530  T=TERRF 

3540  IF  T=1  THEN  CPRINT  "TAPE  ERROR" 

3545  IF  T=1  THEN  8349 
3560  CPRINT  "COMPLETE" 

8565  D=BACKF 
3530  GOSUB  6380 
3609  GOTO  5009 

3000  CPRINT  "THIS  PART  OF  PROGRAM  NOT  COMPLETED" 
3010  GOTO  5140 


Figure  A-3  (Continued) 
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:33; , R<283; 


-£j  SUP  SLICE  FILE-; 


i  '-C. 

.“.—  I  T  .“>■£• 

»J2>L_  X  VC 

rr  .'i  \  _ 
x  r  n  /  — 


3  233 

n*  4  u  Q  1 


«-  ♦  <la  t 

Ji r  a; HI /  * 

THEN  LET 
THEN  LET 

Hi;-i 


R(X>=ft''1338 


V 

;  ■•vv 


4  r  o  ✓ 
r  r 
x  r 

r  p 

IF  fl+Al»-2 


then  let  Cs  rt ; =3i/  1333 
THEM  LET  M=* 


THEM  213 


2o3 

-»A 
X.  »  V 


PRIMT  “CORRECT I MG  BASEL I ME  DRIFT- 
LET  P=R( 1 ) 
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Figure  A-4  Complete  BASIC  Deconvolution  Program  for  Perkin-Elmer 
Instrument 
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.« .”1  •> 
*»Oi. 

LET  31=31+1 

.«  A  0 
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LET  32=32+1 

4  AA 
*t  70 

PRINT  TflB<fll>;  “*** 

CAA 
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caa 

■JX.U 

NEXT  X 

533 
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543 
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Figure  A-4  (Continued) 
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